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Abstract 

We present a comprehensive overview of an extension of the Standard 
Model that contains three right-handed (sterile) neutrinos with masses below 
the electroweak scale [the Neutrino Minimal Standard Model, (i/MSM)]. We 
consider the history of the Universe from the inflationary era through today 



Q H i and demonstrate that most of the observed phenomena beyond the Standard 



Model can be explained within the framework of this model. We review the 



D | mechanism of baryon asymmetry of the Universe in the z/MSM and discuss 



a dark matter candidate that can be warm or cold and satisfies all existing 
constraints. From the viewpoint of particle physics the model provides an 
explanation for neutrino flavor oscillations. Verification of the z^MSM is 
possible with existing experimental techniques. 
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1 Introduction 

^ . The Standard Model (SM) of elementary particles [H, Q, 0] , defined as a renormal- 

izable field theory, based on the SU(3)xSU(2)xU(l) gauge group, and containing 
three fermionic families — left-handed particles, SU(2) doublets, right-handed par- 
ticles, SU(2) singlets (no right-handed neutrinos) and one Higgs doublet — has 
successfully predicted a number of particles and their properties. However, there is 
no doubt that the SM is not a final theory. Indeed, over the past several decades it 
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has become increasingly clear that it fails to explain a number of observed phenom- 
ena in particle physics, astrophysics, and cosmology. These phenomena beyond the 
SM (BSM) are (a) neutrino oscillations (transition between neutrinos of different 
flavors) , (b ) baryon asymmetry (excess of matter over anti-matter in the Universe) , 
(c) dark matter (about 80% of all matter in the Universe consisting of unknown 
particles), (d) inflation (a period of rapid accelerated expansion in the early Uni- 
verse) , and (e ) dark energy (late-time accelerated expansion of the Universe) This 
list of well-established observational drawbacks of the SM is considered complete 
at present. All the other BSM problems — for, instance, the gauge hierarchy and 
strong-CP problems — require theoretical fine-tuning. 

At what energies should the SM be superseded by some other, more funda- 
mental theory? The existence of gravity with the coupling related to the Planck 
scale — Mpi = G^ 1 ^ 2 = 1.2 x 10 19 GeV, where Gat is the Newtonian gravitational 
constant — implies that the cut-off is at least below the Planck scale. If the cutoff 
is identified with Mpi, the low-energy Lagrangian can contain all sorts of higher- 
dimensional SU(3)xSU(2)xU(l)-invariant operators that are suppressed by the 
Planck scale, 



where Csm is the Lagrangian of the SM. These operators lead to a number of 
physical effects that cannot be described by the SM, such as neutrino masses and 
mixings, proton decay, etc. However, as we will discuss below, even the theory 
shown in Eq. (CQ) does not survive when confronted with different experiments and 
observations. 

Alternatively, one can place a cut-off A <C Mpi in Eq. ([T]), which would im- 
ply that new physics (and new particles) appears well below the Planck scale at 
energies E ~ A. If A ^> Mw, where Mw is the mass of the weak W boson, the 
resulting theory suffers from the so-called gauge hierarchy problem, that is, the 
problem of quantum stability of the mass of the Higgs boson against quantum 
corrections from heavy particles. 

Most of the research in BSM physics carried out during the past few decades 
was devoted to solving the gauge hierarchy problem. Many different suggestions 
were proposed concerning how to achieve the "naturalness" of electroweak sym- 
metry breaking. These propositions are based on supersymmetry, technicolor, and 
large extra dimensions, among other ideas. Finding a solution to the gauge hier- 
archy problem, coupled with the need to solve observational and other fine-tuning 
problems of the SM, is extremely challenging. Most of the approaches postulate 
the existence of new particles with masses above the electroweak scale (ranging 
from 10 2 GeV to 10 15 -10 16 GeV). As a result, the proposed theories contain a 
plethora of (not yet observed) new particles and parameters. 
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In this review, we describe a conceptually different scenario for BSM physics 
and its consequences for astrophysics and cosmology in an attempt to address 
the BSM problems named above without introducing new energy scales (that is, 
in addition to the electroweak and the Planck scales). In such an approach, the 
hierarchy problem is shifted to the Planck scale, and there is no reason to believe 
that the field theoretical logic is still a_p_plicable to it. 

Below we show (following Refs. jl, [5j and a number of subsequent works) that 
this goal may be achieved with a very simple extension of the SM. The only new 
particles, added to the SM Lagrangian are three gauge-singlet fermions (i.e., sterile 
neutrinos) with masses below the electroweak scale. Right-handed neutrinos are 
strongly motivated by the observation of neutrino flavor oscillations. In Section [2] 
we review neutrino oscillations and introduce the corresponding Lagrangian. We 
summarize the choice of parameters of the Neutrino Minimal Standard Model 
(r/MSM) in Section [3J In Section HI we present a z/MSM cosmology. We discuss 
the restrictions from astrophysics, cosmology, and particle physics experiments, as 
well as future searches in Section [5j In Section [BJ we conclude with a discussion of 
possible extensions of the i/MSM and potential astrophysical applications of sterile 
neutrinos. 



2 Neutrino oscillations 

When the SM was conceived, neutrinos were thought to be massless and different 
lepton numbers were believed to be conserved (which would not be the case if 
right-handed neutrinos were present). This was a reason for not introducing right- 
handed neutrinos decades ago. 

Transitions between neutrinos of different flavors have been observed. These 
neutrino oscillations manifest themselves in experiments with solar, atmospheric, 
accelerator and reactor neutrinos. The subject of neutrino oscillations has received 
a great deal of attention in recent years and many reviews are available (see e.g. 
Refs. @, 0). The observed oscillations determine the mass differences between 
propagation states [s[ 

Am s 2 un = 7.65i°js x KT 5 eV 2 (2) 

and 

|Am 2 tm |=2.40^i 2 xl(T 3 eV 2 . (3) 

Importantly, each of these results was confirmed by two different types of experi- 
ments: accelerator [i| for Eq. Q and reactor 10] for Eq. (^). 



The observed pattern of neutrino oscillations cannot be explained by the action 
Eq. (CQ) with the Planck-scale cutoff. Indeed, the lowest-order five-dimensional 
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operator 

O, = A af3 (LJ) (tfL}) (4) 



leads to the Majorana neutrino masses of the order m v ~ v 2 /Mpi ~ eV, where 
L a are left-handed leptonic doublets, the index a = e, fi,T labels generations, 
is a Higgs doublet with = z(t2)*<^, c is the sign of charge conjugation, and 
v = 174 GeV is the vacuum expectation value of the Higgs field. 

The fact that the m v following from this Lagrangian is so small compared with 
the lower bound on neutrino mass arising from the observations of neutrino oscilla- 
tions (m u > a/| Am^ tm | ~ 0.05 eV) rules out the conjecture that the theory shown 
in Eq. (TjQ) is a viable effective field theory up to the Planck scale. Therefore, the 
existence of neutrino oscillations requires adding new particles to the Lagrangian 
of the Standard Model. 

Let us add M right-handed neutrinos Nj (I = 1, A/"). The most general renor- 
malizable Lagrangian has the form: 

C = C SM + iNjd^Nj - Z a Nj4> - ^ NfNj + h.c)j , (5) 

where F a i are new Yukawa couplings. The right-handed neutrinos have zero elec- 
tric, weak and strong charges; therefore, they are often termed "singlet," or "ster- 
ile" fermions. The Majorana masses Mj are consistent with the gauge symmetries 
of the SM. Without loss of generality, the Majorana mass matrix in diagonal form 
can be chosen. 

If the Dirac masses Mo = F n r((p) are much smaller than the Majorana masses 



Mi, the type I seesaw formula [Hi, [12|, [13|, [14j holds that 

(m„U = ~ S(^) Q 7t^-(M£) //3 , (6) 

7=1 1 

where m u is a 3 x 3 matrix of active neutrino masses, mixings, and (possible) 
CP-violating phases. An elementary analysis of Eq. ([6]) shows that the number of 
right-handed singlet fermions M must be at least two to fit the data of neutrino 
oscillations. If there were only one sterile neutrino, then the two active neutrinos 
would be massless. If there were two singlet fermions, only one of the active neu- 
trinos would be massless, which does not contradict the results from experiment. 
Moreover, in this case there are 11 new parameters in the Lagrangian (0) — two 
Majorana masses, two Dirac masses, four mixing angles, and three CP-violating 
phases — which is more than the number of parameters (7) describing the mass 
matrix of active neutrinos with one zero eigenvalue. In other words, for M = 2 the 
Lagrangian (0) can describe the pattern of neutrino masses and mixings observed 
experimentally. Of course, the situation is even more relaxed for M = 3. The 
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Lagrangian (JHJ) with M = 3 restores the symmetry between quarks and leptons: 
Every left quark and lepton has a right counterpart. 

The seesaw formula © allows the mass of singlet neutrinos to be a free pa- 
rameter: Multiplying Mp by any number x and Mj by x 2 does not change the 
right-hand side of the formula. Therefore, the choice of Mj is a matter of theoret- 
ical prejudice that cannot be fixed by active-neutrino experiments alone. 

In this review, our choice of Mj is roughly of the order of the other mass term in 
the Lagrangian of the SM, the mass of the Higgs boson. This choice does not lead 
to any intermediate scale between the electroweak and Planck scales, but it does 
require small Yukawa couplings F a j. It allows us to solve, in a unified manner, the 
observational problems of the SM discussed above. Specifically, the parameters of 
the model (jSJ) can be chosen such that it provides a mechanism to generate baryon 
asymmetry of the Universe (BAU) and a dark matter candidate in the form of 
the lightest sterile neutrino. Also, this theory can accommodate inflation, if one 
considers non-minimal coupling between the Hig gs boson and gravity 15( . Finally, 
a scale-invariant extension of the model la, LL7J that includes unimodular gravity 
may solve the problem of stability of the Higgs mass against radiative corrections, 
even if the Planck scale is included, and could lead to an explanation of dark 
energy and of the absence of the cosmo logical constant. 

The Lagrangian (J5J with such a choice of parameters, explaining all confirmed 
BSM phenomenajl] has been termed the z/MSM — the Neutrino Minimal Standard 
Model 0,13. 



3 A choice of parameters: the z/MSM 

In this section we discuss in detail how the parameters of the Lagrangian (jSJ) should 
be chosen so as to explain the BSM phenomena without introducing new physics 
above the electroweak scale. 



3.1 Neutrino masses and oscillations 

As discussed above, the Lagrangian (jHJ) can explain any pattern of active neutrino 
masses and mixings. A naturalness argument 23J, reproduced below, indicates 



■"■There are also several unexplained phenomena in astrophysics (e.g., the 0.511-MeV line 
from the Galactic center [13], pulsar- kick velocities [l9j], a feature in the positron spectrum 
at ~ 100 GeV [2(1 , etc.). We do not consider them in this review, as for these phenomena 
standard physics explanations may exist. Note, however, that the sterile neutrinos may play an 
important role in some of these cases (see Section [5]). In addition, several anomalies in particle 



physics experiments (such as evidence for the neutrinoless double-beta decay [21( and the annual 
modulation of signal, observed by DAMA/LIBRA [22|) have been reported. However, none of 
them has yet been confirmed by other experiments. 
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that the masses Mi should be smaller than 10 7 GeV. Indeed, the stability of the 
Higgs mass against radiative corrections, which arise from ^/-exchange in the 
loops, requires 

al 

where Mn is the Higgs mass. Combining this equation with the seesaw formula 
(H leads to 




~ 10 7 GeV . (8) 

Any Mj mass below this value is "natural". Also note that setting Mi to zero 
increases the symmetry of the Lagrangian (introducing lepton number conserva- 
tion) and thus this choice is also stable against radiative corrections. In the z/MSM 
approach, we set Mi below the electroweak scale, so that these new particles are 
on the same footing as known quarks and leptons. 

The observed existence of two mass splittings (solar and atmospheric, Eqs. 
E])) requires at least two sterile neutrinos. To estimate what neutrino data imply 
for the Yukawa couplings, we take the larger of the two mass splittings |Amj[ tm | 
and find from the seesaw relations (H that 



|2 „ VlAmLJM, 9v1f1 _ l5 M, 

~ ^ ~ z X 1U — — 



^ V_J atm^ _ 2 x 10 _ 15 _^ (Q) 

v z GeV 

where \F\ 2 is a typical value of Yukawa couplings F a j. The condition Mi < 
10 2 GeV implies that |F| 2 < 1(T 13 . Clearly, the theory given in Eq. (jSJ) can 
describe neutrino oscillation data even if one of the three sterile neutrinos has 
arbitrarily small Yukawa couplings. 

Another possible choice of parameters of the Lagrangian ©, suffering from 



the hierarchy problem described in Reference [23|], relates the smallness of the 
neutrino masses to the largeness of the Majorana mass terms Mi. Indeed, if 
we take the Yukawa couplings F a i ~ 1 and the Majorana masses in the range 
Mi ~ 10 10 — 10 15 GeV, we obtain the neutrino masses, as required by experi- 
mental observations. An attractive feature of this scenario is that this new scale 
may be associated with the Grand Unification Theory (GUT) scale. A model (j5 ) 
with this choice of parameters can also give rise to BAU through leptogenesis J24 1 



and anomalous electroweak number non-conservation at high temperatures 
However, we do not discuss this possibility here; for review of the GUT-scale see- 
saw and the thermal leptogenesis scenario associated with it, see, for instance, 
Reference (261 @ 



2 There is still another choice of seesaw parameters with a sterile neutrino mass ~ 1 cV [271 ] . 
which we do not discuss here. 
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3.2 Dark matter candidate 



It has been noticed that a sterile neutrino may make an interesting dark matter 
candidate 28|, |29j, |30j, |3l|, |32j. In the i^MSM, a dark matter sterile neutrino is 
simply one of the singlet fermions (for definiteness, we consider it to be Ni). The 
interaction strength between the sterile neutrino and the matter is superweak with 
the characteristic strength 9Gf, where Gf is the Fermi constant and where the 
mixing angle 1 is defined as 



i' 2 IF, 



al 



a=e,fi,T 



M{ 



(10) 



Based on the universal Tremaine-Gunn bound 33] , the masses of sterile neutrinos 
are restricted to the keV range and above. Specifically, for fermionic DM parti- 
cles an average phase-space density in any dark matter-dominated system cannot 
exceed the density given by the Pauli exclusion principle. Applied to the small- 
est dark matter-dominated objects [dwarf spheroidal galaxies of the Milky Way, 



34 



(dSphs)] this bound translates into M DM > 400 eV 

A sterile neutrino is an example of decaying DM. Through its mixing with the 
ordinary neutrinos, N\ can decay (via Z boson exchange) into three (anti)neutri- 
nos. To be DM, the lifetime of N\ should be greater than the age of the Universe, 
which restricts the mixing angle. A significantly stronger constraint comes from a 
subdominant one-loop decay channel into a photon and an active neutrino. The 
energy of the produced photon is E 7 = Mi/ '2, and the decay width is given by [35 
36 



9 a G 2 F 
1024tt 4 



sin 2 (2^ 1 )M 1 5 ~ 5.5 x 10~ 22 6{ 



Mi 
keV 



'IT 



where a is the fine-structure constant. The expression for the width of the dom- 
inant decay channel has the same parametric form, with the numeric coefficient 
being ~ 128 times larger. This radiative decay would produce a narrow line, for 
instance, in the diffuse X/7-ray background 
the mixing angle, 



30l . l32| . This gives a restriction on 



Q\ < 1.8 x 10~ 5 ( 



keV 



(12) 



meaning that the lifetime of sterile neutrino DM tn 1 ->2, u exceeds 10 24 sec [3j 
six orders of magnitude longer than the age of the Universe. We discuss X-ray 
constraints in detail in Section 15.11 

The contribution of the N\ sterile neutrino to the masses of the active neutrinos 
is of the order Sm u ~ 9\M 1 . The lifetime constraints mean that for M\ > 2 keV the 
correction 5m u is less than the error bar in the solar neutrino mass difference 
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Figure 1: Constraints on the masses and mixing angles of the dark matter sterile 
neutrino N\ (a) and of two heavier sterile neutrinos N 2j 3 (b). These constraints 
come from astrophysics, cosmology, and neutrino oscillation experiments. Abbre- 
viations: BAU, baryon asymmetry of the Universe; BBN, big bang nucleosynthesis. 



This implies that the neutrino N\ cannot contribute significantly to the active 
neutrino mass matrix. As we show below this conclusion is valid for all admissible 
masses and mixing angles. 

Thus, at least three sterile neutrinos are required to account for DM and explain 
neutrino oscillations. In spite of the very weak coupling of DM sterile neutrino 
N\ with matter, it could have been produced in sufficient amounts in the early 
Universe (see Section POl below) . 



3.3 Two heavier neutrinos 

Two other sterile neutrinos (N 2 ,Ns) should interact with the Standard Model 
particles more strongly than Ni to explain the observed pattern of neutrino oscil- 
lations. Here we outline properties of N 2 and N 3 upon which we elaborate further 
in Section HI 

The masses of these two sterile neutrinos should lie in the range 150 MeV < 
M 2 ^ < 100 GeV and should be degenerate (AM 2j 3 -C M 2 ^). These characteristics 
are required for baryon asymmetry generation (Section 14.21) . In order not to affect 
the predictions of big bang nucleosynthesis, we choose their masses and Yukawa 
couplings so that their lifetime r < 0.1 s (i.e. they decay at temperatures T > 
3 MeV) . These bounds are summarized in Figure HJ 

The presence of these sterile neutrinos can affect DM production. If these 
particles decay before the DM is produced, they can generate significant lepton 
asymmetry. In the presence of lepton asymmetry sterile neutrino production pro- 
ceeds differently 29]. If N 2 3 decay after DM production, they heat the primordial 



plasma, increasing its entropy by a factor of 5* and changing its temperature by 
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S 1 / 3 > 1. As the DM sterile neutrinos decouple from the primeval plasma, they 
take no part in this entropy release, making their average momentum smaller than 
that of the active neutrinos by the same factor, S 1 ^ 3 . As a result, the concentra- 
tion of DM particles, compared to that of active neutrinos, decreases by a factor 
of S, and therefore the mixing angle, required to produce the same amount of DM 
particles today, becomes S times bigger for the same DM mass (compared to the 
S = 1 case) [38|] . Therefore, the computation of sterile neutrino DM production 
and relation between the mass and mixing angle should be performed in the full 
vMSM , taking into account the influence of the neutrinos N 2 ,N 3 . 



4 The z/MSM cosmology 

In this section we describe the history of the Universe in the z/MSM from inflation 
onward. 



4.1 Inflation in the z/MSM 



It has been shown 15| that it is possible to incorporate inflation within the i/MSM 
approach, that is without introducing new particles above the electroweak scale. 
The role of the inflaton is played by the SM Higgs field 0. The non-minimal 
coupling of to the Ricci scalar R 



(13) 



leads to a sufficiently flat effective potential for large values of <fi. The constant £ is 
fixed by the Higgs mass and by the amplitude of scalar fluctuations, obtained from 
COBE's observations of the cosmic microwave background (CMB). For the SM 
model to be a consistent field theory all the way up to the Planck scale, the mass 



of the Higgs boson must lie in the interval 126 GeV < M H < 194 GeV [15j, |39_|, |40 



(see also [411. 142]). For uncertainties in these estimates see [40J]. After inflation the 
Universe heats up to the temperature T = T rc h > 1.5 x 10 13 GeV, creating all of 
the particles included in the SM. However, the singlet fermions of the i/MSM are 
not created either during inflation or durin g re heating of the Universe because of 



the small values of their Yukawa couplings [43 



4.2 Baryon asymmetry of the Universe 

Let us consider the stage of the Universe's evolution between the reheating tem- 
perature of ~ 10 13 GeV and the temperatures of ~ 100 GeV. Our Universe is 
baryon asymmetric - it does not contain antimatter in amounts comparable with 
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matter, to wit: 



VB = V* = (8.8 ± 0.2) x 10- 11 (14) 



where s is the entropy density [441 ] . To produce baryon asymmetry of the Universe 



three conditions should be satisfied 45[: (a) Conservation of the baryon number 



should be violated, (b) CP symmetry should be broken, and (c) the correspond- 
ing processes should be out of thermal equilibrium. Although the SM fulfills all 
there requirements [25|, it cannot produce BAU because there is no first-order 



electroweak phase transition with experimentally allowed Higgs boson masses [46 
In addition, the CP violation in Cabibbo-Kobayashi-Maskawa mixing of quarks 
is unlikely to lead to the observed value of BAU ( TT4"1 ) . 



The only source of baryon number non-conservation in the vMSM (as in the 



SM) is the electroweak anomaly [25|. The field configurations important for baryon 



number non-conservation are known as sphalerons [47j. The sum of baryon and 
lepton numbers is not conserved at temperatures above T sph ~ 100 GeV [25j, and 
the BAU can be generated from the lepton asymmetry at T > T sp ^. 

A detailed description of the system of singlet leptons and active fermions in 
the early Universe is necessarily quite complicated. The number of relevant zero- 
temperature degrees of freedom (three active neutrinos, three sterile neutrinos, and 
their antiparticles) is large, and the timescales of different processes can vary by 
many orders of magnitude. To understand the qualitative picture of the creation 
of the matter excess (see Refs. 0, HgJ for a detailed quantitative discussion and 49 



for the original proposal of baryogenesis in singlet-fermion oscillations), one must 
determine the reasons for departure from thermal equilibrium j45|. The singlet 
fermions, with their very weak couplings to the SM fields, play an important role. 
The production rate of iVj is of the order of 

rj ~ |F 7 | 2 T , (15) 

where |F/| 2 = J2 a \Fai\ 2 - Because at T re h singlet fermions were absent, the number 
of created singlet fermions at temperature T is of the order of 

q{T) = = — , is 

n eq T 

where n eq is an equilibrium concentration at given temperature. The temperature- 
time relation is given by t = M ~ Mpi/1.66^/g e ff, where g e ff is the number 
of effectively massless degrees of freedom. The relation shown in Eq. (Tl6l) is valid 
only if q < 1. The singlet fermion equilibrates at temperature T + , defined via 
q(T + ) ~ 1, as 

T + ~ |F 7 | 2 M . (17) 
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The phase-space density constraints, together with X-ray bounds (Section I3.2j) . 
place an upper bound on the Yukawa coupling of the N% at \Fi\ 2 < 10 -21 . The 
relation ffTUj) shows that Ni is irrelevant for baryogenesis because <C 1 at all 
temperatures T > T sp h. In other words, N\ practically decouples from the plasma. 

If for both N 2 and g(T sp h) > 1, N 2 and A3 come to thermal equilibrium 
before T sp h. In this case, even if baryon asymmetry of the Universe were generated 
before equilibration, it would be destroyed in equilibrium reactions with partici- 
pation of sphalerons and singlet fermions. If <7(T sp h) ~ 1, the fermions approach 
thermal equilibrium at the electroweak epoch, exactly what is needed for effective 
baryogenesis. This requirement restricts the masses of A 2i 3 to the 1 — 100-GeV 
region. A detailed analysis of this condition, which accounts for numerical factors 



and complicated flavor structure of asymmetry has been performed 48J . It leads 
to a constraint (Figure [T|) on the mass-mixing angle plane. 

Let us turn now to the effects of CP violation. Although the Yukawa matrix of 
singlet fermions generally contains CP-violating phases, in the tree-level processes 
their effects cancel out. Any quantum loop, which is necessary for CP violation to 
show up, is parametrically suppressed by F 2a , F 3 2 Q , which are much smaller than 
t]b (cf. Eq. ([9])). The only means of overcoming this factor is to have a resonance, 
which occurs when the masses of sterile neutrinos A 2 and A 3 are nearly degenerate. 

Let AM(T) be the mass difference between A 2 and A 3 (temperature dependent 
due to Higgs' vev T dependence and radiative corrections at high temperatures). 
The frequency of oscillations between A 2 and A3 is given by 



10 



\Ml-Ml\ M 2 AM(T) 



E! T 



where we take into account that the typical energy of a singlet fermion Ej ~ T 
and that AM(T) <C M 2 m M3. If the rate of oscillations exceeds the rate of the 
Universe's expansion H(T), the system is out of resonance, and no amplification 
of CP violation occurs. In the opposite case, CP violation has no time to develop. 
This allows us to determine the temperature Tb at which baryogenesis occurs: 
lu ~ H(T B ), which leads to 0, fl, 



T B ~ (M 7 AM(T)M o y /3 . (19) 

The described mechanism is quite effective. The maximal baryon asymmetry A# = 
1 (i.e. 7] B ~ 1CT 2 ) is produced when T B ~ T sph ~ T + , leading to 



n B +n s 



AM(Tb) ~ 0.01 eV ^ 1 V ) • ft ^ s important that interactions with plasma do not 
destroy the quantum- mechanical nature of oscillations of A r 2 3 , in other words that 



< T 



B- 
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The overall qualitative picture is as follows: Coherent pairs of singlet fermions 
are constantly created in interactions with the electroweak plasma. These fermions 
oscillate in a CP violating manner. As a result the Universe develops lepton 
asymmetry in the active, left-handed sector. This asymmetry is converted to 
baryon asymmetry due to anomalous electroweak reactions with baryon number 
non- conservat ion . 

Probably, the simplest way to deal with all these effects in a quantitative way 
is to use the equation for the density matrix p 51, 49, In this case p is a 
12 x 12 matrix (12 = 3x2x2 degrees of freedom for all active and sterile neutrino 
states). For the problem at hand, this equation can be simplified (for details see 
Refs. EQ): 

1. The rates of interactions between active neutrinos are much higher than the rate 
of the Universe's expansion. Therefore, coherent effects for active neutrinos are 
not essential, and the part of the general density matrix p related to active 
leptonic flavors can be replaced with equilibrium concentrations characterized 
by three dimensionless chemical potentials p a (the ordinary chemical potential 
divided by the temperature), giving the leptonic asymmetr y in each flavor. 

2. Active neutrinos acquire temperature dependent masses [52J that are much 
larger than those of singlet fermions. Therefore, all non-diagonal elements of 
the density matrix simultaneously involving the active and sterile states can be 
set to zero. 

3. The DM neutrino Ni is completely decoupled from the system. 

This leaves us with the 2x2 density matrix p^ for singlet fermions iV2 and N$, 
the charge conjugated density matrix p^ for corresponding antiparticles (opposite- 
chirality states, to be precise), and three chemical potentials p a . The correspond- 
ing equations can be written as [H, Hs| : 



.dp 



N 



dt 

.dpN 
>,— — 
dt 

i—r~ 
dt 



[H, p N ] - - {T N , p N - p eq ] + ip a T% , 

[H* 7 p N ] ~ 1 ~{T* N , p N - p eq } - ip a f%* , and 



-ir a L p a + iTr f a L (p N -p ec <) 



— iTr 



(20) 
(21) 
(22) 



Here, Ti, = p(t) + 7i + 7i int is the Hermitian effective Hamiltonian incorporating 
the medium effects on neutrino propagation; p(t) is the neutrino momentum, with 
(p{t)) ~ 3T (we assume that all the neutral fermion masses are much smaller 



than the temperature); Hq 



M 2 
2p(t) 



(we include AMjj to Ti. 



intj] P 



"1 



exp (-p/T) 



is an equilibrium diagonal density matrix; and [,]({, }) corresponds to the 
commutator (anti-commutator). In Eq. (|22|) there is no summation over a and 



1 Li 



via 



are real. The explicit expressions for the equilibration rates T^, T^, 
Yukawa couplings can be found in Ref. 0] for the case in which the temperature 
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is higher than the electroweak scale and in Ref. [48( for the case in which the 
temperature is smaller. These equilibration rates are all related to the absorptive 
parts of the two point functions for active or sterile neutrino states, and they 
contain a square of Yukawa couplings F a j. The real parts of the corresponding 
graphs, together with the mass-squared difference between A2 and A3 determine 
the effective Hamiltonian 7i. For high temperatures T > T sp h the equilibration 
processes are associated with Higgs, W and Z decays to singlet and active fermions, 
to corresponding inverse processes, and to tt — > Nv scattering (where t is the top- 
quark). At smaller temperatures T < T sp h the rates are associated with W and 
Z exchange and with singlet-active mixing through the Higgs vev. The last terms 
in Eqs. ( 1201) and pTj), as well as Eq. ( 1221) are crucial for the generation of lepton 



asymmetry (see p, |48|). In the earlier work (49| the computations and qualitative 
discussion were based on incomplete kinetic equations, that lacked these terms. 

Equations (J20H221), supplemented by initial conditions p^ = Pn = Ha = 
fixed by inflation (Section 14.11) . give a basis for the analysis of the baryon and 



lepton asymmetry generation. It has been shown in Refs. p, |48J that the observed 
baryon asymmetry ( Tl4l) can be generated in a wide range of parameters of the 
i/MSM, provided iV^ are sufficiently degenerate: 

M2 \ J 800, normal hierarchy 



v ' V 1 dtml \GeV J \ 6400, inverse hierarchy v ' 

The generation of baryon asymmetry stops at the temperature of the sphaleron 
freeze-out T sp h- The generation of the lepton asymmetry may still occur at T < 
Tsph [48j |- Moreover, the magnitude of low-temperature lepton asymmetry has 
nothing to do with baryon asymmetry and may exceed it by many orders of mag- 
nitude. Sufficiently large lepton asymmetry = ^ e ~ n ^ e > 10~ 3 can affect the 
production of the DM sterile neutrino (see Section 

Two different mechanisms can lead to late lepton asymmetry generation. As 
discussed above, the singlet fermions enter into thermal equilibrium at T + . As 
the temperature decreases below the electroweak scale, the rate of interactions of 
singlet fermions with plasma gets suppressed by the Fermi constant, and these 
fermions decouple at some temperature T_, that is somewhat higher than their 
mass 48j. The lepton asymmetry can then be generated below T_ in their oscil- 
lations, just as the baryon asymmetry was created above the electroweak scale. 
Alternatively, the lepton asymmetry can be created via decays of A^. The anal- 
ysis ji^] shows that large lepton asymmetries (A^ > 10 -3 ) can be produced, pro- 
vided that the zero-temperature mass degeneracy of A^ is high enough: AM (0) < 
10" 4 v / ]Arr4J (note, that AM(0) can differ significantly from AM(T B ).) The 

3 We denote by n„ e (n,p c ) the number density of the electron neutrino (antineutrino) , and 
assume throughout the paper that the lepton asymmetry is flavor blind. 
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maximal possible lepton asymmetry that can be produced at T_ or during the 
decays of A^ 2 ,3 is A™ ax = 4/(9 x 2 + 4) = 2/11, where 4 is the total number of 
spin-states of iV 2j 3 and 9 is the number of spin-states of three lepton generations. 
This degree of asymmetry can be attained if CP violation is maximal. Generation 
of the baryon asymmetry at the level shown in Eq. ( Tl4l) . together with the large 
lepton asymmetry, is possible for the singlet fermion masses in the 0(GeV) range 
and \Fo\ 



|i*3|. Following the Refs. |48l. l53| we characterize the lepton asymmetry 
as L 6 = 10 6 (n^ - n Pe )/s. The value A™ ax = 2/11 corresponds to L™ x « 700. 
This value is still considerably smaller than the upper limit, imposed by the big 



bang nucleosynthesis [54( L 



2500. 



4.3 DM production in the z/MSM 



In the i^MSM DM sterile neutrinos are produced in the early Universe due to 
their coupling to active neutrinos. An estimate of the rate of DM sterile neutrino 



production Tjv at temperatures below the electroweak scale is given by [55 

T N ~T u e 2 M {T) , 



(24) 



where T u ~ G 2 F T 5 is the rate of active neutrino production, and where 9m(T) is a 
temperature- (and momentum-) dependent mixing angle: 

o\ 



o{ 



e 2 



(T) 



l+*(b(p,T)±c(T) 



(25) 



0\ 



Here 56 



b(p, T) = ^p(2 + cos 2 M^Sr, <T) = 3V2G F (l + sin 2 6 W ) [n h 
naw o60 \ / 



7vr 2 T 4 



n L 



iraw 



(26) 

where 9\y is the weak mixing angle, aw is the weak coupling constant, and p ~ 
few T is the typical momentum of created DM neutrinos. The term c(T) in 025p 
is proportional to the lepton asymmetry and contributes with the opposite sign to 
the mixing of Ni with active neutrinos and antineutrinos. 

If the term b(p,T) dominates c(T) for p ~ (2 — 3)T [which we refer to as 
non- resonant production NRP), the production rate (12~4"]) is strongly suppressed at 
temperatures above a few hundred MeV and peaks roughly at [28| 



T, 



peak 



130 



Mi 
1 keV 



1/3 



MeV , 



(27) 



corresponding to the temperature of the quantum chromo dynamics crossover for 
keV scale sterile neutrinos. This complicates the analysis, because neither the 
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quark, nor the hadron description of plasma is applicable at these temperatures 57 ]. 
In the region of the parameter space (#1, Mi), as demonstrated by X-ray observa- 
tions together with dSph constraints (Section [3. 2p T N (T peak ) <C H(T peak ). In other 
words the DM sterile neutrinos were never in thermal equilibrium in the early Uni- 
verse. This fact allows us to make first principle computation of their abundance 
via certain equilibrium finite temperature Green's functions computed in the SM 



without inclusion of right-handed fermions [57|, |58|, |53j. (See Refs. [3jJ, |59j, [60 
for earlier computations of sterile neutrino DM abundance. The sterile neutrino 
oscillations in the medium were also considered in Ref. [6l|). 

The shape of the primordial momentum distribution of the NRP sterile neu- 
trinos is roughly proportional to that of the active neutrinos (30l |: 



fm(t,p) 



X 



?p/T v {t) 



(28) 



where the normalization x ~ $1 an d where T u (t) is the temperature of the active 
neutrinos. The true distribution of sterile neutrinos is in fact colder than that 
shown in Eq. (128]) . Specifically, the maximum of p 2 /iVi(p) occurs at p/T u pa 1.5 — 1.8 
depending on Mi), as compared with p pa 2.2T U for the case shown in Eq. 



57, 58 



The NRP mechanism specifies a minimal amount of DM that will be produced 
for given Ml and 9\. The requirement that 100% of DM be produced via NRP 
places an upper bound on the mixing angle 9\ for a given mass. This conclusion can 
only be affected by entropy dilution arising from the decay of some heavy particles 
below the temperatures given in Eq. (|2*7|) 381 ] . In the i^MSM such dilution may 



occur due to the late decay of two heavier sterile neutrinos (see Section l3~3l . 

The production of sterile neutrino DM may substantially change in the presence 
of lepton asymmetry. If the denominator in Eq. ( 1251) is small, 



1p 



b(p,T)±c(T) 
Ml 



0. 



(29) 



then resonant production (RP) of sterile neutrinos 29j occurs, analogous to the 



Mikheyev-Smirnov-Wolfenstein effect [62j, [63j. In this case the dispersion relations 



for active and sterile neutrinos cross each other at some momentum p, determined 
from Eq. ( 1291) . This results in the effective transfer of an excess of active neutrinos 
(or antineutrinos) to the population of DM sterile neutrinos. The maximal amount 
of sterile neutrino DM that can be produced in such a way (n res ) is limited by 
the value of lepton asymmetry: n res < \n Ue — n„ e \. The present DM abundance 
DM ~ 0.2 |44j translates into n res /s < 10~ 6 ^, in other words the resonant effects 
in production of sterile neutrinos are not essential for L$ < 1 (for masses in the 
keV range). 
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Resonantly produced DM sterile neutrinos have spectra that significantly differ 



from those in the NRP case [29|, |53j. As a rule, the primordial velocity distribution 
of RP sterile neutrinos contains narrow resonant (cold) and nonresonant (warm) 
components. The spectra for Mi = 3 keV for L 6 = 10, 16, 25 are shown in Fig. [2k 
as functions of comoving momentum q = p/T u . The non-resonant tails are self- 
similar for the same mass and for different lepton asymmetries, and for q > 3 they 



behave as rescaled NRP spectra [58| with the same mass. We term this rescaling 
coefficient the fraction of the NRP component / NRP . The dependence of the average 
momentum (q) on the mass and lepton asymmetry is shown in Fig. [2b. 

4.4 Big Bang Nucleosynthesis. 

Decays of heavy singlet fermions A^ 2 ,3 heat the primeval plasma and distorts the 
spectra of active neutrinos from their thermal values. If this distortion occurs 
around the time of BBN, it may significantly influence its predictions. It has been 



shown [6J, [65J that if the lifetime of sterile neutrinos T2,3 < 0.1 s (i.e. they decay at 
temperatures above ~ 3 MeV), the decay products have enough time to thermalize 
without spoiling the nucleosynthesis. This requirement places a lower bound on 
the strength of Yukawa interactions F a2 and F a3 (Fig. [T]) [66( . A detailed treatment 
of the decays of sterile neutrinos and their interaction with primordial plasma can 
somewhat change these constraints, in particular by allowing for smaller Yukawa 
couplings. 

4.5 Structure formation in the z/MSM 



Comparing the production temperatures Eq. I jzTj) of DM sterile neutrinos with 
their masses shows that they are produced relativistically in the radiation-dominated 
epoch. Relativistic particles stream out of the overdense regions and erase pri- 
mordial density fluctuations at scales below the free- streaming horizon (FSH) - 
particles' horizon when they becomes nonrelativistic. (For a detailed discussion 



of characteristic scales see e.g. 67J] and references therein). This effect influ- 
ences the formation of structures. If DM particles decouple nonrelativistically 
(cold DM models, CDM) the structure formation occurs in a "bottom-up" man- 
ner: specifically, smaller scale objects form first and then merge into the larger 
ones [68J. CDM models fit modern cosmological data well. In the case of parti- 
cles, produced relativistically and remaining relativistic into the matter-dominated 
epoch (i.e. hot DM, HDM), the structure formation goes in a "top-down" fash- 
ion [s^j, where the first structures to collapse have sizes comparable to the Hub- 
ble size 0, H Q- 

The HDM scenarios contradict large-scale structure (LSS) 



observations. Sterile neutrino DM are produced out of equilibrium at the tem- 
peratures shown in Eq. (}2"Tj) and is then redshifted to nonrelativistic velocities in 
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Figure 2: (a) Characteristic form of the resonantly produce (RP) sterile neutrino 
spectra for Mi = 3 keV and L§ = 16 (red solid line) and its resonant and non- 
resonant components. Also shown are spectra for L$ = 10 and Lq = 25, all of 
which have the same form for q > 3. All spectra are normalized to give the 
same DM abundance. For the spectrum with L 6 = 16, its nonresonantly pro- 
duce (NRP) component / NRP ~ 0.12 ((red dashed- dotted line) and its colder, RP 
component (red dashed line) are shown separately. The spectrum with Lq = 10 
has f NRP ~ 0.53 and L 6 = 25 has f NRP ~ 0.60. The RP components happen to 
peak around g res ~ 0.25 — 1. Their width and height depends on L 6 and M x . (b) 
Dependence of the ratio of the average momentum of sterile neutrinos to that of 
the active neutrinos on mass M x and lepton asymmetry. The black solid curve 
represents the ratio of average momenta, computed over the exact NRP spectra to 
the ratio using the approximate analytic form shown in Eq. ( 1281) . The black solid 
curve is the ratio of average momenta, computed over exact NRP spectra to the 
ratio, computed using the approximate analytic form (1281) . For each mass there 
is lepton asymmetry for which the (q) reaches its minimum (q) m i n ~ 0.3{q) Ua . 
The corresponding spectra are the most distinct from the NRP spectra. The frac- 
tion / NRP for the former can be < 20%. 
permission. 



Panel b reproduced from Ref. 53| with 
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the radiation-dominated epoch (warm DM, WDM). Structure formation in WDM 
models is similar to that in CDM models at distances above the free streaming 
scale. Below this scale density fluctuations are suppressed, compared with the 



CDM case. The free-streaming scale can be estimated as [71 



where 1 Mpc is the (comoving) horizon at the time when momentum of active 
neutrinos (p u ) ~ 1 keV. If the spectrum of sterile neutrinos is nonthermal, then the 
non-relativistic transition and A F ° shift by (pn) / (Pu)- To account for the influence 
of primordial velocities on the evolution of density perturbations, it is convenient 
to introduce a transfer function (TF) with respect to the CDM model, 

T(k) EE [Pwdm(£0/Pcom(£0] 1/2 < 1 • (31) 

where the power spectra -P W dm(&) an d -Pcdm(^) are computed with the same cos- 
mological parameters. In general, the characteristic shape of a TF depends upon 
a particular WDM model. In the z/MSM the TF's shape is determined by the 
mass of DM particle and lepton asymmetry, which define the primordial velocity 
distribution. In the case of both RP and NRP sterile neutrino, the TF starts to 
deviated from one for k > k® SH , where 



^FSH 7-T-T7 I [ ^ I 



1 keV J V h 



0.1 \ / Mi 



n Af h 2 J V 1 keV 



h/Mpc . 

(32) 

is a free- streaming horizon for a model with mass Mi and the velocity distribution 
in Eq. (EHD (c.f. 0). In Eq. ( 1321) Qm is the total matter density, and H = 



h x 100 km s 1 Mpc 1 is the Hubble constant today. For k ^> k® m the TF of 
NRP sterile neutrinos falls off as ~ A; -10 (c.f. (z3, S Ifij). For the TF of RP 
sterile neutrinos we define k^ ~ k® m (q)/q rcs . For scales k < k^ the T RP (k) is 
indistinguishable from that of cold-plus-warm dark matter (the mixture of the 
WDM component in the form of NRP sterile neutrinos and pure CDM) with the 
same mass and the same fraction of the warm component F WDM = / NRP . For 
smaller scales (k > k F f H ) the T RP (k) decays, because the cold component of RP 
sterile neutrinos also has a nonnegligible free-streaming, determined roughly by 
fcpsg. This characteristic behavior of NRP and RP spectra is demonstrated in 
Fig. [3j We discuss cosmological restrictions on the sterile neutrino DM models in 
Section 15.1.31 
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Figure 3: Transfer functions (TFs) for the nonresonant production (NRP) spectra 
where M\ = 14 keV and Mi = 3 keV. Also shown is a resonant production (RP) 
spectrum where Mi = 3 keV and L 6 = 16 together with a cold/warm DM spectrum 
for M = 3 keV and 

iwDM — 0.2 (blue dashed-dotted line). Vertical lines mark the 
free-streaming horizon of the nonresonant (fcp SH , left line) and resonant [k r F f H , right 
line) components. Abbreviations: CWDM, cold-plus-warm dark matter; WDM, 
warm dark matter. 
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5 Search for z/MSM physics 

5.1 Combined Restrictions from Astrophysics 

Below we summarize the astrophysical bounds on the sterile neutrino dark matter. 

5.1.1 Phase-space density bounds. 

The most robust lower bound on the DM mass comes from the analysis of the 
phase-space density of compact objects such as dwarf spheroidal galaxies of the 



Milky Way. The universal Tremaine-Gunn bound [33j can be strengthened if one 
(a) knows the primordial velocity distribution of the DM particles and (b) takes 
into account that the dynamics of DM particles is collisionless and dissipationless 
and that an initial phase-space volume occupied by DM particles is therefore pre- 



served in the course of evolution. (For further discussion, see, e.g., |34l. l75l. \7&\ and 
references therein). For the NRP sterile neutrino, this leads to M\ > 1.8 keV 34 . 
For the same mass, the primordial velocity distribution of RP sterile neutrino DM 
is colder than the NRP dark matter. Analysis of available spectra shows that mod- 



els in which M\ > 1 keV are allowed for most values of lepton asymmetries [34 
Note, however, that in derivation of these bounds (unlike the bound, based on the 
Pauli exclusion principle) the presence of dissipative baryons was not taken into 
account. However, because dSphs are the most dark and DM-dominated objects 



(with a mass-to-light ratio > iOOM Q /L ), we expect the results of Ref. [34j to be 
robust. 



5.1.2 X-ray constraints. 



Because the sterile neutrino is a decaying DM candidate, its decays would produce 



a narrow line in spectra of DM-dominated astrophysical objects 30|, I32J. The 



width of this line is determined by the virial motion of DM particles in halos. It 
ranges from AE/E ~ 10 -2 for galaxy clusters to 10~ 4 for dSphs. The expected 
flux of the DM decay is given by 



Mgr e, 

AnDj Mi 



6.38 



M 



fov 



10 10 M ( 







Mpc 



x sin 2 (2^ 



Mi 
keV 



1 5 



keV 



cm' 



sec 



(33) 

where M^J is the mass of DM within a telescope's field of view (FoV), the decay 
rate Y is given by Eq. fl33|) and D L is the luminous distance to the object. Eq. fl33l is 
valid in nearly all interesting cases apart from that of the signal from the Milky Way 
halo (c.f. 77J, |78j). We begin our discussion of X-ray constraints with considering 
the mass range 0.4 keV < Mi < 30 keV with corresponding photon energy E 1 = 
Mi/2 falling into the X-ray band. Modern X-ray instruments have a FoV of 
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~ 10' — 15' which is below the characteristic scale of DM distribution in a large 
variety of objects. As a result, the DM density does not change significantly 
within the instrument's FoV, and the DM flux f[3"3l is determined by the DM 
column density S = J p DU (r)dr: 



fov 



57T 

line of sight 



pDu(r)dr . 



(34) 



Note, that in this limit the expected signal does not depend on the distance to the 
object. Instead, it is determined by the angular size f2f ov and the DM overdensity 
S. The central DM column density is comparable (within a factor ~ 10) for 
most of DM-dominated objects of different types and scales. The objects with 
the highest column densities are central regions of dSphs (such as Ursa Minor or 
Draco) and of Andromeda galaxy. The expected signals from these objects are a 
few times stronger than that of the Milky Way halo (through which one observes 
all the objects). However there are a variety of objects, covering the FoV of an 
X-ray instrument, that would produce comparable signals (in contrast with signals 
from annihilating DM). This gives us a certain freedom to choose observational 
targets, that may be exploited in two ways (described below). Galaxies and galaxy 
clusters are X-ray bright objects whose spectra contain many atomic lines. The 
DM decay line should be distinguished from these lines (and from those of the 
instrumental origin). One can reliably identify atomic lines in cases where several 
of them originate from transitions in the same element and therefore have known 
relative intensities. The emission of hot gas (including the positions and intensities 
of lines) can be described by models with a few parameters such as abundances of 
elements. Another way to distinguish an unidentified astrophysical line from the 
DM decay line is to study its surface-brightness profile. The surface brightness 
profiles of atomic lines are proportional to the gas density squared, while the 
surface brightness profile of a DM line is proportional to the density. To avoid 
complicated modeling and related uncertainties, one can analyze dark outskirts 
of extended objects (such as clusters 79[ or big galaxies [80j) or search for a 
line against a featureless astrophysical background such as the extragalactic X-ray 



background |8jJ, |73, |78|, |82[ . The best strategy is to look in the spectra of X-ray 
quiet and DM dominated objects, in particular dSphs 73, |78| . Extensive searches 
for these lines have been performed using XMM-Newton 18 ll 

INTEGRAL ' 



80 , Chandra 85, 83, 86 



82 




83|, [781, 184J 



87, 88 , and Suzaku 



They did 



not reveal any "candidate" line in the energy range from ~ 0.5 keV to about 
14 MeV. These searches yielded an upper bound on the possible DM flux. To 
convert it to the upper bound on the mixing angle (via Eq. (TTTj) ) , one must reliably 
determine the amount of DM in the observed objects. Determining the amount 
of DM in any given object is subject to specific systematic uncertainties, related 
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e.g. to assumptions about the spherical symmetry of the DM halo, modeling of 
the contribution of baryons to the rotation curves, and assumptions about the 
hydrostatic equilibrium of intra-cluster gas, etc (see e.g. 



SalSfl, 1231). Therefore, it 



is important to rule out the same parameter values from several targets (preferably 
of different scale and nature). 

The bounds obtained from the combined results of these searches (Figure H|) 
are based on independent results from various objects (extragalactic X-ray back- 



ground |8lLl77|, the Milky Way [73, ]85 



78 
89 



Si 182, lS7|, ISSl, M31 [SJ, ISDJ, galaxy 



clusters [79j, |86j], and dSph galaxies [7Sl, |89j). Superposing several results makes 
them robust. Comparing the of production curves (i.e. the relations between 8\ 
and Mi, which produce the required DM abundance) with the bounds on DM de- 
cay flux leads to the upper mass bound. In the NRP case, the upper bound on the 
DM sterile neutrino mass was found to be Mi < 4 keV SQ|. For a more effective 
RP scenario, the smaller mixing angles are admissible, and the corresponding mass 
bound is much weaker: Mi < 50 keV 53j, |34J . 

To improve the modern astrophysical bounds (Fig.HJ) with existing instruments, 
one should conduct prolonged (several Msec) observations of dSph galaxies and 
other X-ray quiet objects, then combine various observations. The bounds would 
improve roughly as Texp with cumulative exposure time. 

A real improvement in sensitivity can be achieved with the next generation of 
X-ray spectrometers which will have an energy resolution AE/E ~ 10 -3 -10~ 4 and 
which will maximize the product of their effective area and FoV. For a comparison 
of the sensitivities of existing and future missions, see Ref. 0,0] . For an example 
of a spectrometer with optimal characteristics for performing a decaying DM search 
see e.g. 92|, |9l|. Such a mission would allow us to improve the sensitivity towards 
the DM search by up to a factor of 10 3 . Its preferred observational targets would 
be our Galaxy, its dwarf satellites, and galaxy clusters. 

If a candidate line is found, its surface-brightness profile can be measured, dis- 
tinguished from astrophysical lines, and compared with those of several objects 
with the same expected signal. This would make an astrophysical search for de- 
caying DM another type of a direct detection experiment. 



5.1.3 Lyman-a constraints. 

Below we describe the Lyman-alpha method and the Lyman-alpha constraints on 
the properties of warm and cold-plus-warm dark matter models and in particular 
sterile neutrino dark matter. 



Structure formation at different scales. The sterile neutrino DM creates a 
cut-off in the power spectrum of matter density fluctuations at sub-Mpc scales. 
At larger scales, the power spectrum can be determined by CMB experiments, 
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Figure 4: The allowed region of parameters for DM sterile neutrinos produced 
via mixing with the active neutrinos (unshaded region). The two thick black 
lines bounding this region represent production curves for nonresonant produc- 
tion (NRP) (upper line, L 6 = 0) and for resonant production (RP) (lower line, 



T max 



u 



700) with the maximal lepton asymmetry, attainable in the z/MSM [5 
The thin colored curves between these lines represent production curves for (from 
top to bottom) Lq = 8,12,16,25, and 70. The red shaded region in the upper 



right corner represents X-ray constraints 77], |78j, |80j, |88|, [89] (rescaled by a factor 
of two to account for possible systematic uncertainties in the determination of 
DM content 



86 



80] ). The black dashed-dotted line approximately shows the RP 
models with minimal (q) for each mass, i.e., the family of models with the largest 
cold component. The black filled circles along this line are compatible with the 
Lyman-a bounds 9(|, and the points with M\ < 4 keV are also compatible with 
X-ray bounds. The region below 1 keV is ruled out according to the phase-space 
density arguments 34J. Abbreviation: BBN, big bang nucleosynthesis. 
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or reconstructed from observed galaxy power spectrum (via LSS surveys). For 
any mass and mixing angle that satisfies the constraints named above, the sterile 
neutrino DM model fits the CMB and LSS data as well as the ACDM "concordance 
model" . One of the main probes of the matter power spectrum at sub-Mpc scales 
is the Lyman-a forest method, discussed below. 



Lyman-a method. It has been well established by analytical calculations and 
hydrodynamical simulations that Lyman-a absorption lines in spectra of distant 
quasars are produced by clouds of neutral hydrogen at different redshifts along the 
line of sight. This neutral hydrogen is part of the warm (~ 10 4 -K) and photoionized 
intergalactic medium. Opacity fluctuations in the spectra arise from fluctuations in 
the neutral hydrogen density, which can be used to trace cosmological fluctuations 
on scales k ~ (0.1 — 5) h Mpc" 1 , at redshifts z ~ 2 — 4. Thus it is possible to 



94 . l95| . For each quasar, 



infer fluctuations in the total matter distribution [! 
the observed spectrum I(z) can be expanded in (one-dimensional) Fourier space. 
The expectation value of the squared Fourier spectrum is known as the flux power 
spectrum Pp(fc). 

The flux power spectrum Pf(/c) is a complicated function of the cosmological 
parameters. In the range of scales probed by Lyman-a data, the effects of free- 
streaming may be compensated for by a change in other cosmological parameters. 
Therefore, it is important to perform collective fits to the Lyman-a data and to 
CMB and LSS data. At redshifts and scales probed by Lyman-a, the fluctuations 
already enter into a mildly non-linear stage of gravitational collapse (Sp/p > 1). 
To predict Pf(^) for a given cosmological model, one must perform hydrodynami- 
cal simulations: for CMB experiments probing the linear matter power spectrum, 
it is sufficient to compute the evolution of linear cosmological perturbations using 
a Boltzmann code such as camb [96j]). Hydrodynamical simulations are necessary 
both for simulating the non-linear stage of structure formation and for computing 
the evolution of thermodynamical quantities (so as relate the non-linear matter 
power spectrum to the observable flux power spectrum). In principle, to fit Lyman- 
a data one should perform a full hydrodynamical simulation for each cosmological 
model. This is computationally prohibitive, and various s implifying; approxima- 



tions have been proposed [97j, |98|, |99|, [100|, UM, UM, UM l!Q3|. Although potentially 



very powerful, the Lyman-a method is indirect and hinges on a large number of 
assumptions (see e.g. [671] ). 



Lyman-a bounds. A Lyman-a analysis 104L \74 . 1 1051 . 1 1061 . 11071 ] has been per- 



formed for NRP mechanism of st erile neut rinos, assuming a model with only one 
sterile neutrino. The bounds 74j, Il06. Il05lj were recently revised in Ref. |67lj. who 

and 



used the same SDSS Lyman-a data 



99 combined with WMAP5 data 
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paying special attention to the interpretation of statistics in the parameter extrac- 
tion and to possible systematic uncertainties. It was shown that a conservative 
(i.e. frequentist, 3-er) lower bound is Mi > 8 keV. The Lyman-a method is still 
under development, and there is a possibility that some of the related physical 
processes are not yet fully understood. However, it is difficult to identify a source 
of uncertainty that could give rise to systematic errors affecting the result by more 
than 30%. Even with such an uncertainty, the possibility of all DM being in the 
form of NRP sterile neutrinos can be ruled out by comparing Lyman-a results 
with X-ray upper bounds [67 . 

A full Lyman-a analysis of DM models containing both resonant and non- 
resonant components has not yet been conducted. However in 67j a cold-plus- 
warm DM model (CWDM) has been analyzed. Although the phase-space distri- 
bution of RP sterile neutrinos does not coincide exactly with such mixed models, 
some results can be inferred from the CWDM analysis. In particular, it has been 



shown [90j, that the cosmological signature of RP sterile neutrino DM can be ap- 
proximated by that of CWDM. With these results 67], it was shown that for each 
mass Mi > 2 keV, there exists at least one model of sterile neutrino that can 
account for the totality of dark matter, and that is consistent with Lyman-a and 
other cosmological data. These corresponding values of lepton asymmetry can be 
obtained within the z/MSM. 



5.1.4 Galaxy substructures. 



The suppression of the matter power spectrum at sub-Mpc scales may affect galaxy 
substructures. Qualitatively one expects a cut-off scale for the existence of DM 
substructures in the Milky Way-type galaxies. The primordial velocities of DM 
particles do not allow them to move too close together, thereby creating shal- 
low density profiles. Quantitative conclusions about the mass function of dwarf 
satellites and slopes of density profiles of galactic substructures require numerical 
simulations, as the galaxy formation today is at a non-linear stage. 

A number of studies have claimed discrepancies between astrophysi cal obser 



108, 109 



yations and CDM numerical simulations of galaxy formation (see e.g 
lid . 111 . 73]). For example: 
1. The number of observed dSphs of the Milky Way is still more than an order 
of magnitude below the CDM predictio ns, in spit e of the drastically improved 



sensitivity in the dSphs search (see e.g. [1121 . Ill3| ): 



2. direct measurement s of DM density profiles in dSphs seem to favor cored profiles 



114 Ilia . 11161 . 1112 



3. there ma y be i ndic ations of the existence of the smallest scale at which DM is 
observed 1121 . Ill7l . 

However, it is not clear whether a "CDM substructure crisis" exists. Various 
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astrophy s ical r esolutions of such a crisis has been proposed (see e.g. [118l . lll9l . ll20 



12ll . Il22l . I123L Il24l . Il25| ). It is extremely challenging to verify these resolutions 



both experimentally and numerically because the full-scale numerical simulations 
including baryons and their feedback on structure formation are not currently 
possible. To understand how galaxy structure formation occurs in the WDM 
models, one must perform numeric al si mulations, that take into account thermal 



velocities of DM particles (see e.g. [1261 ] . [671] ). 



5.2 Particle physics experiments 

Because many of the parameters of the z/MSM are fixed or constrained by exper- 
iments on neutrino masses and oscillations, as well as by cosmological considera- 
tions, it is possible to rule out or to confirm this model. 

We list below the z/MSM predictions for experiments that are under way or 
that will be carried out in the near future. 

1. Nothing but the Higgs boso n w i thin the mass interval 126 GeV < Mh < 
194 GeV will be found at the LHC " 
are given in 



4Q|). 



1271 . l43l . l39l . l40| (uncertainties in these numbers 



2. One of the active neutrinos in the z/MSM must be very light, m 1 < C»(10" 6 ) eV 



|ill28 



which fixes the masses of the two other active neutrinos at m 2 ~ 9 • 10 
eV and 7723 ~ 5 ■ 10~ 2 eV for normal hierarchy or 7712,3 — 5 ■ 10~ 2 eV for the inverted 
hierarchy. As a result, a n eff ective Majorana mass for neutrinoless double-beta de- 
cay can be determined 1291 ] . For normal (inverted) hierarchy the constraints are 
1.3 meV < mpp < 3.4 meV (13 meV < mpp < 50 meV). 

3. There will be negative results for the direct and indirect searches for weakly 
interacting massive particles, proton decay, neutron oscillations, and extra sources 
of CP-violation in the hadronic sector. 

Although the experiments listed above can rule out the z/MSM, they cannot pro- 
vide its direct confirmation. We can only obtain confirmation if all three singlet 
fermions are discovered. Astrophysical searches for the DM sterile neutrino are de- 
scribed at length above. Laboratory searches for this particle would require a chal - 



130 



lenging, detailed analysis of the kinematics of /5-decays of different isotopes 
A pair of new neutral leptons, creating the baryon asymmetry of the universe, 
can be searched for 66j in dedicated experiments that use existing intensive pro- 
ton beams at CERN, FNAL, and planned neutrino facilities in Japan (J- PARC). 
Alternatively, one can look for a specific missing energy signal in K, D, and B 
decays. In general, the z/MSM paradigm tells us that to uncover new phenomena 
in particle physics we should utilize high-intensity proton beams with fixed-target 
experiments or to very high intensity charm or B-factories, rather than high energy 
accelerators. 
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6 Conclusion 

Sterile neutrinos may play a key role in resolving the BSM problems. For the 
specific choice of parameters described above, they can explain the neutrino flavor 
oscillations and the origin of the matter-antimatter asymmetry in the Universe, 
and can provide a DM candidate. Sterile neutrino DM can be cold or warm, 
affecting the formation of structures at sub-Mpc scales. This parameter choice 
does not contradict existing experimental bounds and it assumes masses below 
the electroweak scale. This means that sterile neutrinos can be searched for in 
existing accelerator and laboratory experiments. The parameter space of the sterile 
neutrinos in the z/MSM is bounded from all sides, allowing for verification or 
falsification of this model. 



Beyond the z/MSM. The z/MSM is a minimal theory that can address the SM 
problems in a unified way. Of course, the fact that it is minimal does not mean 
that it is correct. Even if the low-energy electroweak theory indeed contains three 
relatively light neutral right-handed leptonic states, the full theory may be much 
more complicated. Extending the z/MSM at low or high energies will change the 
predictions of this model. For example, introducing the low-energy supersymmetry 
would add another DM candidate to the theory, allowing wider variation of the 
parameters of DM sterile neutrinos. If GUT theories are proven to be true, they 
would introduce another source of baryon number non-conservation, leading to a 
different origin of BAU. 

Finally, we briefly discuss two modest modifications of the z/MSM. In the first, 
one adds to the Lagrangian (jHJ) the higher-dimensional operators given in Eq. (Tj[|). 
This allows for production of DM sterile neutrinos at reheating after inflation 



431 ]. provided that they are heavier than 100 keV. This does not contradict to 
X-ray constraints because the mixing angle Q\ is not essential for this production 
mechanism. The sterile neutrino, in this case, can only play the role of CDM 
candidate. 

T he se cond modification of the z/MSM requires the addition of a singlet scalar 



field jX 3 II ] . This field may give a mass to singlet fermions due t o Yu kawa couplings 



to these fermions, and it may play the role of the inflaton 13l|, |l32|. The decays of 



this scalar to DM sterile neutrinos may serve as an efficient mecha nism of cold or 



warm DM production, without conflicting with X-ray constraints [13ll . Il33l . [134 



135. 136 . In such models the mass of sterile neutrino DM can be in the MeV 



range or even in the GeV range. The radiative decays of such DM particles can be 
searched for with Fermi (cf. jl37| ). whose large FoV requires a new search strategy. 
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Astrophysical applications of sterile neutrinos. A final word about astro- 
physical applications of right-handed neutrinos, discussed in this work. The mean 
free path of a sterile neutrino inside a collapsing star is much larger than that 
of an ordinary neutrino. The escape of sterile neutrinos from the supernova may 
thus provide an effective cooling mechanism. For sterile neutrinos with the masses 
M > 10 2 keV the analysis of the data from the supernova 1987 A leads to an uppe r 



bound on the mixing with electron neutrino sin 2 (2# e ) < 10 10 1381 . Il39l . Il40l . Il41 
In the case of sterile neutrinos with the mass M < 10 2 keV the analysis shows 
that the mixing with the ele ctron neu trinos is bounded by sin 2 (26 l e ) < f-f^-) 
for 1 keV < M < 100 keVjpl [HI, 

while the mixing with v^ T is bounded by 
sin 2 (2^ )T ) < 3 x 10- 8 T 3 o 6 H^c.fnSlI) where T 30 = T core /30 MeV and a typical 
supernovae core temperatures, T core , lie in the range 20 — 70 MeV [1431 ] . Owing to 
large uncertainties in the supernova physics, the constraints on the mixing angle 
of the sterile neutrino that arise from these considerations happen to be weaker 



65 



than the X-ray bounds discussed abo ve (S e ction I57L21) [64 . 

At the same time, the analysis of 14J, 1 1451 ] shows that sterile neutrinos with 
the masses 1 keV < M < 5 keV and with the mixing angle with v e 10~ 10 < 
sin 2 (26* e ) < 10~ 8 may enhance lepton number, energy, and entropy transport in 
supernovae explosions. The results of 1461 ] indicate that sterile neutrinos with the 
mass M 2j3 ~ 200 MeV and mixing angles 6> 2 > 10~ 8 or 10 -8 < # 2 < 10~ 7 can 
augment shock energies of core collapse supernovae. The asymmetric emission of 
keV sterile neutrinos from a supernova explos i on m ay le ad t o an explanation of 
pulsar kick velocities 147 . 148 . 149 . 133 . 134 . 150j |. see 151 ] for a review. The 
X-ray photons from early decays of these particles can influence the formation of 
molecular hydrogen and can be im portant for th e early star fo rmation and re- 



ionization 



152, 153, 154, 155, 156 



157. 158. 159 for a discussion 



See also Refs. 

of the formation of super-massive black holes and degenerate heavy neutrino stars. 

Finally, if DM consists of many components (as may occur in extensions of 
the z/MSM) such that sterile neutrinos constitute only a fraction / < 1 of it, the 
X-ray bounds shown in Fig. 0] become we aker by the same factor. Moreover, the 
Lyman-a constraints also become diluted 160l . 1671 ] , if the rest of DM is cold. The 



analysis in 67J shows that the DM sterile neutrino with the spectrum given by 
Eq. (1281) and a mass Mi > 5 keV for / < 0.6 satisfies all existing constraints at 
99.7% CL. Once / < 1, sterile neutrino effects in astrophysics may be even more 
important. 
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